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ABSTRACT 

The Resonant Switch (RS) model of twin high-frequency quasi-periodic oscillations (HF QPOs) 
observed in neutron star binary systems, based on switch of the twin oscillations at a resonant 
point, has been applied to the atoll source 4U 1636-53 under assumption that the neutron star 
exterior can be approximated by the Kerr geometry. Strong restrictions of the neutron star pa¬ 
rameters M (mass) and a (spin) arise due to fitting the frequency pairs admitted by the RS model 
to the observed data in the regions related to the resonant points. The most precise variants of 
the RS model are those combining the relativistic precession frequency relations with their mod¬ 
ifications. Here, the neutron star mass and spin estimates given by the RS model are confronted 
with a variety of equations of state (EoS) governing structure of neutron stars in the framework of 
the Hartle-Thorne theory of rotating neutron stars applied for the observationally given rotation 
frequency f^t ~ 580 Hz (or alternatively f^t ~ 290 Hz) of the neutron star at 4U 1636-53. It is 
shown that only two variants of the RS model based on the Kerr approximation are compatible 
with two EoS applied in the Hartle-Thorne theory for /rot ~ 580 Hz, while no variant of the RS 
model is compatible for /ot ~ 290 Hz. The two compatible variants of the RS model are those 
giving the best fits of the observational data. However, a self-consistency test by fitting the ob¬ 
servational data to the RS model with oscillation frequencies governed by the Hartle-Thorne 
geometry described by three spacetime parameters M, a and (quadrupole moment) q related by 
the two available EoS puts strong restrictions. The test admits only one variant of the RS model 
of twin HF QPOs for the Hartle-Thorne theory with the lOandolfi et all ( l2010h EoS predicting 
the parameters of the neutron star M ~ 2.10Mq, a ~ 0.208, and qla^ ~ 1.77. 
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1 Introduction 

The high-frequency quasi-periodic oscillations (HF QPOs) in the Galactic Low Mass 
X-Ray Binaries (LMXBs) containing neutron (quark) stars are often demonstrated as 
two simultaneously observed pairs of peaks (twin peaks) in the Fourier power spec¬ 
tra corresponding to oscillations at the upper and lower frequencies (vu, vl) that sub¬ 
stantially change over time (even in one observational sequence). Most of the twin 
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HF QPOs in the so-called atoll sources ( van der Klisll2006 ) have been detected at lower 
frequencies 600 - 800 Hz vs. upper frequencies 900 - 1200 Hz, demonstrating a clus¬ 


tering of the twin HF QPOs frequency ratio around 3:2 dAbratnpwiczejal 


Belloni et al. 2(ffl7 : Torok et al, 2008albl : Montero & Zanotti 2012 : Wang et ^ 
Stefanovll2014l) 


2005 


2013 


For some atoll neutron star sources the upper and lower HF QPO frequencies can 
be traced along the whole observed range, but the probability to detect both QPOs 
simultaneously increases when th e frequency ratio is close to ratio of small natural 
numbers, namely 3:2, 4:3, 5:4 (lT6rokll2009h . The analysis of root-mean-squared- 
amplitude evolution in the group of six atoll sources (4U 1636-53, 4U 1608-52, 
4U 0614-H09,4U 1728-34,4U 1820-30,4U 1735-44) shows that the upper and lower 
HF QPO amplitudes equal each other and alter their dominance while passing r ational 
frequ ency ratios (3 : 2 or 4 : 3) corresponding to the datapoints clustering (ITorok 
20091) . Such an “energy switch effect” can b e well explained in the framework of non¬ 
linear resonant orbital models as shown in Horak et"^ ( 2009h . Another interesting 
phenomenon related to energy of th e twin HF QPOs has been recently demonstrated in 
Mukheriee & Bhattacharyyal ( 120121) . Further, analysis of the twin peak HF QPO ampli¬ 
tudes in the atoll sources (4U 1636-53,4U 1608-52, and 4U 1820-30,4U 1735-44) 
indicates a cut-off at resonant radii corresponding to the frequency ratios 5 :4 and 4: 3 
respectiyely, implying a possibility that the accretion disc inner edge is located at the in- 
nermost resonant ra dius rather than at the innermost stable circular geodesic (ISCO, 
Stuchlrk et al.l201 1 ). The situation is different for some of the Z-sources where the twin 
peak frequency ratios are clustered close to 2: 1, and 3 :1 ratios as demonstrated in 
the case of Circinus X-1 (IBoutloukos et al.l 120061) . Then the resonant radii could be 


expected at larger distance from the ISCO than in the atoll sources ( Torok et al ]m^. 

The eyolution of the lower and upper twin HF QPOs frequencies in the atoll and Z 
sources suggests a rough agreement of the data distribution with the so-called hot spot 
models of HF QPOs, e.g., the Relatiyistic Precession (RP) model prescribing the eyo¬ 
lution of the upper frequency by the Keplerian freq uency vu - vk and the lowe r fre¬ 
quency by the precession frequency Vl = Vk - 


dStella & Vielrilll998l Il999h 


goy- 

erned by the radial epicyclic frequency of geodetical circular motion. In rough agree¬ 
ment with the data are other models based on the assumption of the oscillatory mo¬ 
tion of hot spots, or accretion disc oscillations, with oscillatory frequencies giyen by 
the geodetica l orbital and epicyclic motion. Th ey include the modi hed RPl model 
( Bursal 2005 ). the To tal Precession (TP ) model ( Stuchlrk et al.l 2007 ), the Tidal Dis- 
ruption (TD ) model (IKostic et al.l 120091) . or the Warp Disc oscillations (WD) model 
( KatolEo^ ). In all of them the frequency difference vu - vl decreases with increasing 
magnitude of the lower and upper frequencies, in accord with the obseryational data 
( Belloni et alJl2007 : Barret et d]|2005at Torok et d]|2012 ). This property of the obser- 
yational data excludes the epicyclic oscillations model assuming vu = vg and vl = 
(Urbanecet^ 2010bl) that works well in the case of HF QPOs in black hole LMXBs 
jTordk et al.ll2005hPI 


^In the black hole sources, twin peaks with fixed pair of frequencies at the ratio 3 :2 are usually observed 
and can be explained by the in ternal non-linear resonance of oscillations with radial and vertical epicyclic 
frequencies iTdrok et i3l^. 

^Note tFiat quite recently a special frequency set of HF QPOs has been reported for the neutron star binary 
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The vu/vl frequency relations, given by the models mentioned above, can be com¬ 
pared to the observa tional data found fo r neutron star LMX Bs, e.g., data of the atoll 
source 4U 1636-53 ( Barret et alJ 2005a : Toroket al. 2008alf9l . or the Z-source Circi- 
nus X-1 ( Boutloukos et al.ll2006 ). The parameters of the neutron star spacetime can be 
then determined due to the fits of the data to the frequency-relation mo dels. The rotat¬ 
ing ne utron stars are described properly by the Hartle-Thorne geometry (IHartle & Thorne 
19681) characterized by three parameters: mass M, internal angular momentum J and 
quadrupole moment Q, or by dimensionless parameters a - JjM^ (spin) and q - 
QMjfi. In the special case when qja^ - 1, the Hartle-Thorne external geometry 
reduces to the well known Kerr geometry if it is expanded up to the second order 
in a. The Kerr approximation is very convenient for calculations in strong gravity 
regime because of simplicity of relevant formulae. It has been shown recently that 
near-maximum-mass neutron (quark) star Hartle-Thorne models, constructed for any 
given equation of state (EoS), imply qja^ ~ 1 and the Kerr geometry is applicable with 


high precision in such si tuations instead of the Hartle-Thorne geometry (lUrbanec et al 


20131: Torok et d]|2010h . Such high-mass neutron stars can be expected at the LMXB 
systems due to the mass increase caused by the accretion process. 

Assuming the geodesic orbital and epicyclic frequencies related to the Kerr ge¬ 
ometry, the fitting procedure applied to the RP model of the frequency-relation evo¬ 
lution implies mass-spin relation M{a) — Mo |^1 -i- k(a + fl^)j rather than concrete val¬ 
ues of the neutron star pa rameters M and a ; for the Z-source Circinus X-1 there is 
Mo ~ 2.2Mo and k ~ 0.5 ( Torok et al. 2010l) . The same mass-spin relations, but with 


different values of the Schwarzschild (no-rotation) mass Mq ~ 1.8 Mp and the c on- 
stant k ~ 0.75, were obtained for the atoll source 4U 1636-53 ( Torok et al. 2012 ). In 


the case of the models similar to the RP model (RPl, TP), the same M{a) relations were 
found, while for th e models TD and WD, the relations are different - for details see 
Torok et al. ( 2012h . Quality of the fits to the data obtained for individual m odels is very 


poor f or the atoll source 4U 1636-53. This fact is extensively discussed in iTorok et al 


( 20121). Bad fitting of o bservational data with the frequency-relation models was found 
also in lLin et ah ( 2011 ) for the atoll source 4U 1636-53 and the Z-source Sco X-1 for 
some models of the HE QPOs with the frequency relations given by some phenomena 
of non-ge odesic origin (iMiller et al.lll998t IZhang et alJ 120061 : iMukhopadhvavl 120091 : 

^1201 ih . 

The bad fitting of the data distribution in the atoll sources by the frequency-relation 
models of HE QPOs based on the assumption of the geodesic character of the os¬ 
cillatory frequencies invoked attempts to find a correction of a non-geodesic origin 
reflecting some important physical ingredients, as influence of th e magnetic field of 
the neutron star o nto slightly charged innermost parts of the disc ( Bakalu_et al]|2010l: 

Kovaf et ^ 2008h . of thickness of non -slender oscillating tori ( Sttuub_&_Sramko^ _ 

2009), or of oscillating string loop model ( Stuchlrk & Koloa2012 , 20141 : Cremaschini & Stuchlik 
2013b . Such modifications of the frequency-relation models could make the fitting 


system XTE J1701—407 that is one of the least luminous atoll sources with Lx ~ 0.01 Ledd jPawar et all 
l2013h . This frequency set resembles observations of the HF QPOs in the microquasars, i.e., black hole 
binary systems, and it can be explained by the model of sti'ing loop oscillations IStuchhk & Kolo^2Q15l) that 
works quite well also in the case of Galactic microquasars GRS 1915+105, XTE 1550-564, GRO 1655-40 
IStuchlflt & Koloa2014h . 
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procedure better as shown for a simple toy model in Torok et al. ( 2012h . However, 
in all these cases, some additional free parameter has to be introduced along with 
the spacetime parameters of the neutron star. Some relevant modifications can be 
also obtained in the framework of models related t o the braneworld compact objects 
dStuchlfk & Kotrlovall2009HSchee & Stuchi^l2009h . 

Therefore, the Resonant Switch (RS) model of twin peak HF QPOs has been re¬ 
cently proposed modifying the standard orbital frequency-relation models in a way 
that allows to keep the assumption of the relevant frequencies being combinations of 
the geodesic orbital and epicyclic frequencies. No non-ge odesic corrections are nec- 
essary in the RS model, although these are not excluded ( Stuchlrk et al. 20121 2013 ) 
- the RS model considers only the spacetime parameters of the neutron star exterior 
as free para meters. The RS mod el has been applied in the case of the atoll source 
4U 1636-53 ( Stuchlfk et al.ll201^ and tested for this atoll source by fitting the obser¬ 
vational data using the frequency relat ions predicted by th e RS model as acceptable 
due to the neutron star structure theory ( Stuchlfk et alJl2014 ). 

The fitting procedure predicts the mas s and spin paramete rs of the 4U 1636-53 
neutron star with relatively high precision (IStuchlrk et al.ll2014l) . Here we test the fre¬ 
quency relation pairs of the RS model giving the best fits for the corresponding values 
of the mass M and spin a of the neutron star, using variety of equations of state con¬ 
sidered recently in modeling the rotating neutron stars in the framework of the Hartle- 
Thorne theory. Strong limits implied by the Hartle-Thorne models can be obtained 
due to the precise knowledge of the rotation velocity of the 4U 1636-53 neutron star 
( Strohmaver & Markward)i2002 ). We are then able to put strong restrictions on validity 
of the acceptable frequency-relation variants of the RS model. 


Resonant switch model of twin HF QPOs 
in the 4U 1636-53 atoll source 


2.1 The RS model 


We br i efly su mmarize the basic ideas of the RS model - for details see IStuchlrk et al 
(I 20 1211201^ . According to the RS model a switch of twin oscillatory modes creating 
sequences of the lower and upper HF QPOs occurs at a resonant point. Non-linear 
resonant phenomena are able to excite a new oscillatory mode (or two new oscillatory 
modes) and damp one of the previously acting modes (or both the previous modes)lj 
Switching from one pair of the oscillatory modes to some other pair will be relevant up 
to the following resonant point where the sequence of twin HF QPOs ends. 

Here, two resonant points at the disc radii Xout and Xin are assumed (x = r/(GM/c^) 
is the dimensionless radius, expressed in terms of the gravitational radius), with ob¬ 
served frequencies and v“, v“, being in commensurable ratios : 

and p™ - n™ : m”. Observations of the twin H F QPOs in the atoll systems put 
the restrictions v“ > 


and p™ < ( ToroklE009 ). In the region related to the res- 


^Note that the sw itch could occur fo r other reasons, e.g., due to the phenomena related to the magnetic 
field of neutron stars IZhang et al.l2006l) . 
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onant point at Xout, the twin oscillatory modes with the upper (lower) frequency are 
determined by the functions M, a) (v™X.r; M, a)). In the region related to the in¬ 

ner resonant point at Xin different oscillatory modes given by the frequency functions 
Vy (x; M, a) and v”(x; M, a) occur. All the frequency functions are assumed to be com¬ 
binations of the orbital and epicyclic frequencies of the geodesic circular motion in 
the Kerr backgrounds. Such a simplihcation is correct with high accuracy for neutron 
(quark) stars with large masses, close to maximum allowed for a given equation of state 
( Torok et al. 2010l: Urbanec et al. 2013 ). that can be assumed in the known atoll or Z- 
sources because of mass increasing due to the accretion. Of course, for neutron stars 
having mass signihcantly lower than the maximal allowed mass, the Hartle-Thorne ex¬ 
ternal geometry reflecting al so the role of the quadrupole moments of the neutr on star 
has to be taken into account ( Urbanec et al.ll2013 ; Gondek-Rosiriska et al.ll2014 ). 

In the Kerr spacetime, the vertical epicyclic frequency va and the radial epicyclic 
frequ ency Vr take the form (e.g., Perez et al.ll997 ; Stella & Vietrill998 ; Torok & Stuchlrk 
200'5h 

V^r^arVi, ( 1 ) 


Vfl = as v|, 


where the Keplerian (orbital) frequency vk and the dimensionless quantities determin¬ 
ing the epicyclic frequencies are given by the formulae 


Vk = 

1 (GM 

2n 

ag = 

l_i^ 

x3/2 

Ur = 

1-^.- 


, 1/2 


■J 


3a^ 


x“ 


3/2 


2>a^ 




( 2 ) 

(3) 

(4) 


Detail s of the properties of the or bital and epicyclic frequencies can be found in iTordk & Stuchlrk 
( 2005h : Stuchlrk & Schee ( 2012 ). We can see that any linear combination of the orbital 
and epicyclic frequencies depends equally on the mass parameter M, therefore, their 
frequency ratio becomes independent of M. Then the conditions 


Vy'(x; M, a ): v™'(x; M, a) - p° 
v'^{x-,M,a) : v”(x;M,fl) = p” 


(5) 


imply relations for the spin a in terms of the dimensionless radius x and the resonant 
frequency ratio p that can be expressed as a°p\x) and a“(x), or in an inverse form 

out/ 


Xp{a) and x”(a). 


The frequency-relation functions ha ve to meet the resonant frequencies that can be 
determined by the energy switch effect ( Torok 2009 : Stuchlrk et al. 2012h . In the RS 
model applied here, two resonant points and two pairs of the frequency functions are 
assumed. This enables direct determination of the Kerr background parameters as¬ 
sumed to govern the exterior geometry of the neutron (quark) star. At the resonant 
radii the conditions 


v™'(x; M, a), v” = v‘u(x; M, a) 


(6) 
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are satisfied along the functions M™' 

® Ponl ' ' fin ' ' f II ■ 

(quark) star are then given by the condition (IStuchlrk et al.ll2012li 


(a) an d M™ (a). The para meters of the neutron 




(7) 


This condition predicts M and a with precision implied by the error occurring in deter¬ 
mination of the resonant frequencies by the energy switch effect that is rather h igh for 
the observational data obta ined at present state of the observational devices (see lTorok 
2009t Stuchlfk et alJ 2012h . However, the fitting of the observational data by the fre¬ 
quency relations predicted by the RS model improves substantially the precision of 
determination of the neutron star parameters and, simultaneously, restricts the versions 
of the RS model that can be considered as realistic. 

Starting from the results obtained in Stuchlfk et alJ ( 2012 ). in the present paper we 
consi der pairs of the frequency relations g iven by the RP model (IStella & Vietrilll998 


1999h . the TP model ( Stuchlfk et alJl20()7 ). an d their modihcatio ns RPl (lBursall2005 ). 
and TPl, combined also with the TD model ( Kostic et al.ll2009l) . and the WD model 
( Kato|[2008 ). The frequency relations are summarized in Table [T] In the RS model ap¬ 
plied to the source 4U 1636-53 the frequency relations are co mbined, and the switc h 
of their validity occurs at the outer resonant point as described in lStuchlfk et alJ (l2012h . 
For each of the frequency relations under consideration the frequency resonance func- 
tions and the resonan ce conditions determining the resonant radii ;r„:m(a) are given in 
Stuchlfk et al. ( 2012 ). 


Table 1: Frequency relations corresponding to individual HF QPO models. 


Model 


Relations 


RP 

vl 

1 

> 

II 

Vu 

= Vk 

RPl 

vl 

1 

> 

II 

Vu 

= Vg 

TP 

Vl 

-Vg-Vr 

Vu 

= Vg 

TPl 

Vl 

> 

1 

II 

Vu 

- Vk 

TD 

Vl 

= Vk 

Vu 

= Vk + Vr 

WD 

Vl 

= 2 (vk - Vr) 

Vu 

> 

1 

n 

II 


2.2 Application to the source 4U 1636-53 


The mass M and spin a ranges predicted by the RS model with resonant frequencie s 
given by the energy switch effect are very large (see Table 1 in lStuchlfk et al.ll2012h . 
However, the ranges can be strongly restricted by htting the observational data near 
the resonant points by the pairs of the frequency relations corresponding to the twin 
oscillatory modes. We use the data of tw in HF QPOs in the 4U 1636-53 s ource as pre 


sented a nd studied in iTorok et al.l (12012h . analysed in the original papers bv iBarret et al 


(l2005all^ - in this case it is immediately clear what is the extension of the data related 
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to the resonant points with frequency ratio 3 : 2 and 5 ; 4, respectively. In the fitting pro¬ 
cedure, based on the formulae related to the Kerr spacetime, we applied those switched 
twin frequency relations predicted by the RS model tha t are ac ceptable due to the neu¬ 
tron (quark) star structure theory ( Stuchlrk et al.l 201 3, 20141) . All the resulting twin 
frequency relations considered in our testing are presented in Table |2] where the values 
of the mass M and spin a of the neutron star predicted by the RS model and the re¬ 
lated fitting procedure presented in Stuchlrk et al.l ( 2014 ) are explicitly given along with 
the correspo nding errors. In fi tting the observational data, the standard least-squares 
(X^) method dPress et al.ll2007h has been applied. In the space of the lower and upper 
frequencies, vl and vu, the^^-test represents the minimal (squared) distance of the fre¬ 
quency relation curve vu (vl) given by a model of the twin HF QPOs from the observed 
set of data points. There is 


A 




A„ = Min 


f7=l 


1 \ Pisco 

^n,p 


/ p, 


( 8 ) 


where l„^p is the length of a line between the centroid values of the nth measured 
data point [vL(n), v\j{n)] and a point [vL(p(n)), vu{p(n))] belonging to the relevant fre¬ 
quency curve of the model; the points are considered down to the point corresponding 
to the ISCO. The quantity cr„ ^ denote s the length of the part of this line located within 
the error ellipse around the data point dPress et al.ll2007 ). 

The ^^-test has been applied solely for th e RP, TP and TD fre quency relations along 
the whole range of the observational data in iTorok et al.l d2012h . However, the results 
were quite unsatisfactory, giving ~ 350 and ;^^^/dof ~ 16. On the other hand, 
the RS model enables increase of the fit precision by a lmost one order, giving in the best 


cases ~ 55 and^^/dof ~ 2.6 ( Stuchlrk et al.ll2Q14l) 3 


Th e mass M and spin a ranges determined by the fitting procedure in lStuchlrk et al.l 
( 20141) for acceptable combinations of frequency-relation pairs are illustrated in Fig¬ 
ures [T] and |2] and summarized in Table |2] The best fit is obtained for the combination 
of frequency pairs RPl-i-RP. 

The results of the fitting procedure will be further tested by confrontation with de¬ 
tailed Hartle-Thor ne theoretical models ( Hartle & Thornell968 ; Chandrasekhar & Millet 
I974I: Miller 1977 ) describing slowly rotating neutron stars that are constructed under 
the observationally given constraint of the rotation frequency 580 Hz (o r 290 Hz) rele¬ 
vant for the 4U 1636-53 neutron star ( Strohmaver & Markwardl 20021). using the va - 
riety of widely accepted equations of state that were studied in (lUrbanec et al.ll2013h . 
We assume a detailed test of a much more extended family of acceptable equations of 
state in a future paper. 

The results of the RS model have to be related in future to the limits on the 4U 
1636-53 neutron star parameters indicated by other possible observational phenom¬ 
ena. In fact, a preliminary result of simultaneous treatment of the twin peak HF QPOs 
and profiled (X-ray ) spectral lines indicates the neutron star mass to be M ~ 2.4 M0 
dSanna et al.ll2012h that gives an important restriction on the results of the RS model 


^Tlie fitting procedure has been realized in the ranges of M and a predicted by the RS model with data 
given by the energy switch effect, but we convinced ourselves that outside these ranges the fits are worse 
than inside of them. 
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Table 2: The best fits of the ohservational data and the corresponding spin and mass parameters 
of the neutron star located in the 4U 1636-53 source, along with related errors in determining 
spin and mass of the neutron star due to the fitting procedure. The results are taken, as the most 
promising ones, from lStuchlik et aP ( l2014h . The two cases of the neutron star parameters that 
are in agreement with the Hartle-Thome model of the neutron stars are shaded. They are those 
corresponding to the best fits of the observational data of twin HF QPOs. 


Combination of models 

■ 

^ mm 

a 

Aa 

M [Mq] 

AM [Mq] 

RP1(3:2) -H RP(5;4) 

55 

0.27 

0.02 

2.20 

0.04 

TP(3:2) -1- RP(5;4) 

55 

0.52 

0.02 

2.87 

0.06 

RP1(3:2) -H TP1(5:4) 

61 

0.20 

0.01 

2.12 

0.03 

RP1(3:2) -H TP(5;4) 

62 

0.45 

0.03 

2.46 

0.06 

TP(3:2) -H TP1(5:4) 

68 

0.31 

0.02 

2.39 

0.05 

RP(3:2) -H TP 1(5:4) 

72 

0.46 

0.03 

2.81 

0.09 

WD(3:2) -H TD(5:4) 

113 

0.34 

0.08 

2.84 

0.21 


and restricts substantially the variety of allowed combinations of frequency relations 
used in the RS model. However, we clearly need more detailed study of the profiled 
spectral lines based on the precise predictions of the character of the external spacetime 
of the neutron star. 


3 Hartle-Thorne model of rotating neutron stars 


The Hartle-Thorne theory represents a standard approximative method of construct¬ 
ing models of compact stars (neutron stars, quark stars , white dwarfs) within general 
relativ ity, assuming rigid and slow rotation of the stars ( Hartlel[r967 : Hartle & Thorne! 
Il968l) . It is treating deviations away from the spherical symmetry as perturbations 
with terms up to a specified order of the rotational angular velocity Q of the compact 
star. Going up to second order in Q, the theory gives the lowest order expressions for 
the frame dragging at, the moment of inertia I - JIQ., with J denoting the angular mo¬ 
mentum of the star, the shape distortion caused by centrifugal effects, the quadrupole 
moment Q and the change in the gravitational mass due to rotation 6M. Recent results 
indicate that the slow rotation approach is quite correct for all the observed rotating 
neutron stars, even in the case of the fastest observe d pulsar PSR J1748-2446ad with 
rotational frequency ~ 716 Hz (lUrbanec et al.ll2013l) . For very fast rotation only, near 
to that giving centrifugal break-up, we have to solve numerically the full set of the Ein¬ 
stein equations rather tha n using the approximati ve a pproach of Hartle-T horne theory 


see models presented in Bonazzola et al. (1998) and Stergioulasl ( 2003 ). 


For our purpose s, the second-order slow -rota tion Hartle-Thorne approximat ive 
theory developed in Hartle & Thorne ( 1968 ) and Chandrasekhar & Miller ( 1974 ) is 
quite appropriate because of the rotation frequency observed for the neutron star in 
4U 1636-53. Then the Hartle-Thorne geometry describing both the internal and ex- 
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ternal spacetime takes in the geometric units with c = G = 1 the form 

ds 2 = _e 2 ’'o [1 + 2 /!o(r) + 2 /z 2 (r)P 2 ( 6 ')] df^ 

+e 2 ^o 11 + — [2mo(r) + 2 m 2 (r)P 2 m f dr^ 

+r2 [1 + 2 k 2 (r)P 2 ( 6 ')] {d^^ + [d0 - w(r)df]^ sin^ o], (9) 

where vo, do and the coordinates are identical with corresponding spherical non-rotating 
solution, cx){r) is a perturbation of order fi, representing the frame dragging, and ho{r), 
/j 2 (r), mo(r), m 2 (r), k 2 (r) are perturbations of order Q?. All of these perturbations are 
functions of the radial coordinate only. The non-spherical angular dependence is deter¬ 
mined by the second-order Legendre polynomial P 2 iS) = 5 (3 cos^ 6 - 1^. 

All the perturbation functions have to be calculated under appropriate boundary 
conditions at the centre and at the surface of the compact star. The second-order pertur¬ 
bations are labeled with a subscript indicating their multipole order: I - 0 for spherical 
perturbations, / = 2 for the quadrupole perturbations representing the deviation away 
from the spherical symmetry. By matching the internal and external solution at the star 
surface, the external parameters of the compact star as measured by distant observers 
can be calculated: the mass M, angular momentum J and quadrupole moment Q that 
fully characterize the external gravitational field in the slow-rotation approximation, if 
one is retaining only perturbations up to second order. 

To construct the internal solution, the Einstein equations are solved with the source 
term given by the energy momentum of a perfect fluid. Rigid rotation of an axisym- 
metric configuration means that the four velocity has components 


U‘ = [- {g„ - 2Qg,^ + Q}g4,^y\'^, = Qt/'. 


( 10 ) 


The derivation of the equations for t he perturbation quantiti e s together with boundary 
condit i ons has b een gi ven in detail in lHartle & Thornel(ll968l) : IChandrasekhar & Miller 
( I 974 I) : Millet ( I 977 I 1 . We will not repeat th is in the present paper, using the same 
procedures as those presented in iMillen (Il977h . 


4 Equations of state 


The crucial ingredient of the compact star models is the equation of state describing 
properties of matter constituting them. Neutron stars are expected to consist of neu¬ 
trons closely packed in /1-equilibrium with protons, electrons and at high densities also 
muons, hyperons, kaons and possibly other particles. Their central densities corre¬ 
spond to microphysics that is not well understood, therefore, they serve as laborato¬ 
ries of nuclear matter under extreme conditions, giving complementary information to 
those obtained in the collider experiments. 

A wide range of approaches for nucleon-nucleon interactions and their role in mod¬ 
eling o f the structure of neutron stars has been used - see the review bv iLattimer & Prakash 
(l2007h . An alternative to the standard neutron star picture is represented by quark stars 
consisting partially or fully from deconflned quarks. The most radical version of this 
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approach is represented by strange stars consisting entirely f rom deconfined quarks 
( Farhi & Jaffell l984 : Haensel et al. 19861: Colpi & Millei|[l9^ ). It is based on the sug¬ 
gestion of IWittenI (119841) that matter consisting of equal numbers of up, down and 
strange quarks represent the absolute ground state of strongly interacting matter. It is 
important to mention that the strange stars have to be bound together by a combination 
of the strong and gravitational forces, in contrast to neutron stars where only the grav¬ 
ity is responsible for the binding. Here we restrict attention on the equations of state 
governing the neutron stars only. 

We consider a set of neutron-star matter equations of state that are based on various 
approaches - following the recent st udy of the neutron sta r properties related to the be¬ 
havior of the quadrupole moment (lUrbanec et ah 2013 ). We give a brief review of 
these equations of state, details can be found in lUrbanec et al.l (120131) and the original 
literature. 

We choose relatively wide set o f Skyrme parameterizations, whose labels are start¬ 
ing with S; see IStone et al.l (l2003h for details of Skyrme potential and diferences be¬ 
tween various parameterizations. We use two variants of APR model based on the vari- 
ational theory refle cting the three body forces and the relativistic boost corrections 
( Akmal et al. 1998 ). APR corresponds to A18 + 6v + UIX*, while APR2 corresponds 
to A18 -H UIX, where relativistic boost corrections are not included. The UBS equa¬ 
tion of state is based on th e relativistic Dirac-Brueckner-Hartree-Fock mean field the¬ 
ory (lUrbanec et al.ll2010ah and correspond to model originally labeled as H. The non- 
relativistic Brue ckner-Hartree-Fo ck theory is represented by the equation of stat e la¬ 
beled as B BB2 (Baldo et al. 1997 ). We also use GlendNH3 ( Glendennind 1985 ) and 
BalbNlHl ( Balberg & Galll 1 9~ ) equation of state including the hyperons at high den- 
sities. FPS is the very well know EoS and has been used very often in the past 


( Lorenz et al.lll993l) and BPAL12 is very soft EoS giving very low maximum mass 


(lBombacilll995l) . Stiff equations of state were recently constructed in the framework 


of the auxiliary field diffusion Monte Carlo technique and are labeled as Gandolfi 
( Gandolfi et al. 2010l) . Our selection of EoS represents wide range of possible models, 
however some of these do not meet current observations dSteiner et al.l2010i Demorest et al 


2010t Antoniadis et al. 2013 ). Here, we focus on the selection of the acceptable vari¬ 


ants of the RS model using the Hartle-Thorne theory of neutron stars that is properly 
applied to the restricted set of the equations of state. 

It is useful to give the maximal values of the neutron star parameters M and a ob¬ 
tained in the framework of different approaches to the equation of state. It follows from 
general relativistic restrict ions that mass of a neutron star cannot exceed Mmax ~ 3 M© 
(IRhoades & Ruffinilll9741) . The realistic equations of state put limit on the maximal 
mass of neutron stars MmaxN ~ 2.8 Mq ( Postnikov et al. 20 id) - the extre mal maxi¬ 
mum MmaxN ~ 2.8 Mq is predicted by the field theory ( Muller & Serotll99^ . The limit 


some special case (Miither et al. 

1987) or bv some Skyrme models. The variational 

approaches (Akmal & Pandharioandelll997; 

Akmal et al. 1998) and other approaches 

(Urbanec et al.||2010a|) allow for Mn,,vM ~ 2.25 M^. 


On the neutron star dimensionless spin the limit of a < flmaxN = 0.7 has been re¬ 
cently reported on the b ase of numerical , non-approximative methods, independently 
on the equation of state ( Lo & LinI I 2 OI ih . The Hartle-Thorne theory can be well ap- 
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plied up to the spin a ~ 0.4 ( Urbanec et al.ll2013 ). 

In the Hartle-Thorne models of rotating neutron stars the spin of the star is lin¬ 
early related to its rotation frequency. The rotation frequency of the neutron star at 
the atoll source 4U 1636-53 has been o bserved at fmt ~ 580 Hz (or fmt ~ 290 Hz, if 
we observe doubled radiating structure. IStrohmaver & Markwardtj 2002 ). Such a ro¬ 


tation frequency is sufficiently low in comparison with the mass shedding frequency, 
and the Hartle-Thorne theory can be applied quite well, predicting spins lower enough 
than the maximally allowed spin. The theory of neutron star structu re then implies for 
a wide variety of re alistic equations of state the spin in the range ( Stergioul^ 2003 : 
Urbanec et alJl2013l) 

0.1 <a <0.4. (11) 

Of course, the upper part of the allowed spin range corresponds to the rotation fre¬ 
quency /rot ~ 580 Hz, while the lower part corresponds to /rot ~ 290 Hz. The related 
restriction on the neutron star (near-extreme) mass reads 


M < 2.8Mo. 


( 12 ) 


Detailed comments on the precission of the Kerr geometry approximating the Hartle- 
Th orne geometry in depend ence on the spacetime parameters a and qja^ can be found 


Stuchlrk & KolosI ( 120151) - for the HF QPO models the differences induced by 


the Kerr approximation could be smaller than five percent for the dimensionless spin 
a < 0.4. 


Testing the RS model by equations of state applied 
in the Hartle-Thorne model of the neutron star 
in the source 4U 1636-53 


5.1 Selection of the relevant variants of the RS model 

The resulting limits on the mass M and spin a of the 4U 1636-53 neutron star implied 
by the data fitting procedure realized in the framework of the RS model of HF QPOs 
are presented in Table taken from IStuchlrk et al.l (120141) and reflected in Figures [T] 
and |2j the precision of the mass and spin estimates is also reflected in Figures [T] 
and |2] The Hartle-Thorne models are constr ucted for a variety o f acceptable equa¬ 
tions of state discussed above (and studied in Urbanec et al. 20131) for both possible 
rotational frequencies of the 4U 1636-53 neutron star. We use nine parameteriza- 
tions of the Skyrme equation of state (SkT5, SkO’, SkO, SLy4, Gs, SkI2, SkI5, SGI, 
SV) and other nine equations of state (UBS, APR, FPS, BBB2, BPAL12, BalbNlHl, 
GlendNH3, Gandolfi, APR2) that well represent the variety of eqations of state. The re¬ 
sults of the Hartle-Thorne model that are calculated for the equations of state under 
consideration are illustrated in the M-a plane in Figure[T]for both the assumed neutron 
star rotation frequencies /rot ~ 290 Hz and /rot ~ 580 Hz. Clearly, all the mass-spin 
dependencies constructed for /rot ~ 290 Hz can be excluded. In Figure |2] we give de¬ 
tailed picture of fitting the mass and spin range implied by the acceptable variants of 
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Fig. 1: Hartle-Thome models of neutron stars with rotation f requency frot ~ 290 Hz (red) 
and /rot ~ 580 Flz (blue), for a variety of EoS considered in lUrbanec et al.l (l2013h: SkT5, 
SkO’, Sk O, SLy4, Gs, SkI2, S kI5, SGI, SV - Skyrme equation s dStone et alj I 2 OO 3 I) UB S 
equation dUrbanec et al. 2010ah. FP S equation jLorenz et al.l 1993I). AP R 1 Akmal et^ll998h . 
BBB2 equation i Baldo et al.lll997h . BPAL1 2 equation 1 Bombacilll995[l . BalbNl Hl equation 
I Balbere & Gallll997ll , Gl endNH3 equatio n 1 Glendenning 19851) . APR2 equation dAkmal et alJ 
Il998h . Gandolfi equation dGandolfi et al.ll201ol) . Each equation of state applied in the Hartle- 
Thome model predicts for a fixed rotational frequency a sequence of stable states represented 
by a curve in the M-a plane; its final point indicates an instability. Clearly, all predicted values 
of the 4U 1636-53 neutron star spacetime parameters M, a related to /ot ~ 290 Hz are located 
outside the predictions of the Hartle-Thome models. 


the RS model by the a(M) curves constructed for the acceptable EoS with the rotation 
frequency/ot ~ 580 Hz. 

We immediately see that no equation of state allows to construct a Hartle-Thorne 
model that can fit the RS model data, if we assume the rotational frequency of the 4U 
1636-53 neutron star /ot ~ 290 Hz. For the rotational frequency /ot ~ 580 Hz, 
the Hartle-Thorne models give very interesting restrictions that are in significant agree¬ 
ment with results of the fitting the HE QPO data in the framework of the RS model. 
The Hartle-Thorne model based on the Skyrme equation of state SV meets with high 
precision the prediction of the RPl+RP version of the RS model that gives the best fit 
to the twin peak HE QPO data observed in the 4U 1636-53 source for the neutron star 
parameters M ~ 2.20 Mq and a ~ 0.270 The Hartle-Thorne model based on the Gan- 

^The same precision of the fit, namely ~ 55, is obtained for the TP+RP version of the RS model. 
However, in this case the predicted mass and spin, M ~ 2.87 Mq and a ~ 0.52, are outside the values 
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Fig. 2: Hartle-Thome models of neutron stars with rotation frequency /mt ~ 580 Hz are illus¬ 
trated for the a{M) dependencies related to the two best fits of the observational data by the RS 
model. The acceptable sequences of equilibrium configurations predicted by the Hartle-Thorne 
model are constructed for two equations of state (Skyrme SV EoS for RPl+RP variant, and 
Gandolfi EoS for the RPl-l-TPl variant). Notice that the variant of the RS model predicting 
the lower parameters of M = 2. 12Mq and a = 0.20 is also very close to the final state given 
by the sequence of equilibrium states predicted by the Skyrme Gs EoS, indicating possibility of 
an instability of the neutron star in the 4U 1636-53 system in near future. 
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dolfi equation of state meets with high precision the prediction of the RPl+TPl ver¬ 
sion of the RS model that gives the second best ht to the observational data of the twin 
HF QPOs in 4U 1636-53 for a neutron star having parameters M ~ 2.12Mo and 
a ~ 0.20. Notice that the variant of the RS model giving the second best ht is in accord 
with another version of the Skyrme equation of state (Gs) which predicts the neutron 
star with parameters M ~ 2.11 Mq and a ~ 0.20. Such a result demonstrates that 
the 4U 1636-53 neutron star could be in a state very close to an instability, as the neu¬ 
tron star mass and spin indicated by the HF QPO data htting procedure can corre¬ 
spond to the hnal state of the evolution of the neutron stars rotating with the frequency 
/rot ~ 580 Hz and governed by the Skyrme equation of state Gs - see Figure|2 Predic¬ 
tions of all the other variants of the RS model are located in the M-a plane at positions 
that are evidently out of the scope of all the equations of state considered in the present 
paper - we can expect that this is true even for all the variants of the presently known 
equations of state. 


5.2 Parameters and the shape of the neutron star 


Shape of isobaric surfaces P - const, and the shape of the neutron star surface P - 0 
are given by 

r(P = const., 6) = ro(P) H- /o(?‘o) + / 2 (^o)^’ 2 (cos 6), (13) 


where rn is th e spherical coordinate and functions /o and ^2 are given by the relations 
(lMillerlll977h 

(14) 


k(r) = 


and 


Anr^P + m(r) 
r[r - 2 m(r)]p 2 
Anr^P + m{r) 

1 


Piir) = -h 2 (r) - . 


(15) 


(16) 


The calculati ons have been realized using the detailed set of equations presented in 


,Milled (Il977h . Then equatorial and polar radii governing the surface shape of the rotat¬ 
ing neutron star read 

Req = Ro{P) + URo) - i/2(-Ro) , (17) 


Rpol = Ro{P) + ^o(Ro) + ^2(Ro) ■ 


(18) 


The agreement of the Hartle-Thorne neutron star models based on the Skyrme 
and Gandolh equations of state with two best hts of the observational data of twin 
HF QPOs observed in the source 4U 1636-53 enables to predict in detail properties 
of the neutron star in this source. Namely, we are able to hnd along with the two 
known parameters, mass M and spin a, also the radius in the equatorial plane P(0 - 
njl) and along the symmetry axis R{6 - 0) and whole the shape of the neutron star 
surface, the moment of inertia /, and quadrupole moment Q and its dimensionless 


acceptable by the neutron star models. 
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form q - QMjJ^. Then we can calculate also the parameter representing compactness 
of the neutron stars in dependence on the latitude 




m 

2M ■ 


(19) 


We can consider the characteristic values of the compactness parameter related to 
the equatorial plane 9 - njl and the symmetry axis 0 = 0. Here we give for sim¬ 
plicity the compactness parameter related to the basical, spherically symmetric model 
that is a starting point of the Hartle-Thorne models, given by 


Co 


Rj) 

2Mo' 


( 20 ) 


The results of the Hartle-Thorne model calculations for all three equations of state 
giving acceptable agreement with the data predicted by the RS model are presented 
in Table [3] For the Skyrme equation of state SV, and the Gandolh equation of state, 
the neutron star parameter s are given for the ma ss parameter corresponding to the mean 
value of the data fitting dStuchlik et al.l 120141) . In the case of the Skyrme equation 
of state Gs, the mass parameter of the neutron star (M ~ 2.11 Mq) corresponds to 
the maximal value predicted by this equation of state, i.e., it gives the instability point 
of the neutron stars governed by this equation of state. This mass parameter is lower 
than the related mean value given by the data fitting, but it falls into the allowed range 
of the mass parameter. 


Table 3: Parameters of neutron stars predicted by the selected equations of state giving Hartle- 
Thorne models compatible with the fitting of the twin HF QPOs observed in 4U 1636-53. 
The mass is fixed to the mean value given hy the fitting procedure with exception of the Skyrme 
equation of state Gs where it corresponds to the final state indicating an instability. The spin 
predicted hy the Hartle-Thorne model fits the range given hy the error determined by the data 
fitting. All models predict very compact neutron stars with low value of the parameter qlcP'. 


EoS 

M [Mo] 

Req [km] 

a 

q 

qja^ 

RolilMo) 

Gandolh 

2.12 

10.75 

0.205 

0.0723 

1.71 

1.72 

Gs 

2.11 

10.84 

0.201 

0.0676 

1.68 

1.75 

SV 

2.20 

13.41 

0.272 

0.1940 

2.60 

2.07 


6 Self-consistency test by the Hartle-Thorne geometry 

Assuming that the external geometry of the 4U 1636-53 neutron star can be approx¬ 
imated by the Kerr geometry, we have found that the two most precise variants of 
the RS model can be fitted by realistic EoS applied in the Hartle-Thorne model of 
slowly rotating neutron stars. For the RPl-i-RP variant, the Skyrme SV EoS predicts 
M ~ 2.20Mo, a ~ 0.272, qla^ ~ 2.60. For the RPl-tTPl variant, the Gandolh EoS 


15 











Fig. 3: Sequence of the Hartle-Thome spacetime parameters implied by the Gandolfi, Skyrme 
Gs and Skyrme SV EoS for the rotational frequency f^ot = 580 Hz. Left: Spin (a) as a function of 
mass (M). Right: Reduced quadrupole moment (qla^) as a function of mass (M). All sequences 
are constructed for the mass interval 1.2 Mq < M < M^.^. 


predicts M ~ 2.12 M©, a ~ 0.205, qja^ ~ 1.71; for this RS variant, also the Skyrme Gs 
EoS gives acceptable values of M ~ 2.11 M©, a ~ 0.201, qla^ ~ 1.68, however, this 
estimated mass parameter is at the maximum allowed mass for the EoS. 

We have to realize now a self-consistency test of the Kerr geometry approximation 
applied in fitting the observational data. We have to check, if fitting the observational 
data by the ^^-test using the orbital and epicyclic frequencies related to the Hartle- 
Thorne geometry with parameters M, a, q governed by the acceptable EoS gives results 
comparable or better than the htting based on the assumption of the Kerr geometry 
approximation. We realize the self-consistency test in the following three steps. 

Eirst, we characterize the sequence of states given by the Hartle-Thorne model 
for the acceptable EoS and the rotational frequency of the neutron star /rot = 580 Hz 
- see Eigure |3] The free parameter is the mass M and we represent the sequence 
of the states by the functions a(M) and qla^{M). For each considered EoS we give 
the Hartle-Thorne model for M > 1.2 M©, and we follow the sequence of allowed 
neutron star states to the limiting values of the parameters M, a, qja^ corresponding 
to the maximal allowed mass for the considered EoS. For each of the EoS, the closest 
approach of the Hartle-Thorne model to the Kerr geometry occurs for the maximal 
mass allowed by given EoS. We can see that closest approach to the Kerr geometry is 
allowed for the Hartle-Thorne model based on the Gandolh EoS, enabling the lowest 
value of qla^ ~ 1.5 for Mmax ~ 2.24M©. On the other hand, the largest difference 
from the Kerr geometry approximation are expected for the Skyrme SV model with 
the lowest value of qla^ ~ 1.9 - the maximal mass is highest for this EoS, reaching 




2.4 Mp 


Second, we define the orbital and epicyclic radial and vertical frequencies of the qua¬ 
sicircular geo desic motion in the Hartle-Thorne geometry, VK(r; M, a, q), Vr(r, M, a, q), 
vg(r; M, a, q) ( Abramowicz et all2003l: Torok et al.l2008bl) . For completeness, we give 
the explicit expressions for the frequencies in the Appendix. 

Third, using the Hartle-Thorne orbital and epicyclic frequencies, we repeat the least 
square (j^f^-test) fitting procedure for the same sample of the observ ational data in 
the 4U 1636-53 atoll source as those considered in our previous paper (IStuchlrk et al 
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20141) . For the self-consistency test, we study only the two selected variants of the RS 
model, RPI-hRP and RPl-i-TPl, along the sequences of the neutron star parameters 
M, a, q related to the allowed EoS and the rotational frequency of the neutron star0 
The results of the fitting procedure are given in Figure |4] for the RPl-i-RP variant 
of the RS model and the Skyrme SV EoS. Along with the htting based on the Hartle- 
Thorne geometry, we repeat for comparison also the results obtained under the Kerr 
approximation of the neutron star external spacetime. The htting procedure implies 
the best ht;^f^ = 101; in comparison to the best ht based on the Kerr approximation 
(x^ = 55, M ~ 2.20 Mq, a ~ 0.272), the mass parameter is shifted to lower value 
of M ~ 2.1 IMo, and higher value of spin a ~ 0.286. Moreover, at the values of 
mass and spin predicted by the Kerr approximation (when - 55, Stuchlfk et alJ 
the Hartle-Thorne geometry implies x^ > 1000. Such a large discrepancy 

3 when large errors of 


2014 . 


1000 . 

is caused by relatively large value of the parameter qja^ 


the Kerr approximation are expected. The resulting value of the Hartle-Thorne best ht, 
X^ - 101, is too high in comparison with the Kerr approximation value, = 55, and 
we can conclude that the RPl-i-RP variant of the RS model is not satisfying the self- 
consistency test. 

For the RPl-i-TPl variant of the RS model and the Gandolh EoS, the results of 
the htting procedure are presented in Figure|5] We again give for comparison the results 
of the best ht obtained for the Kerr geometry approximation (x^ - 61, M ~ 2.12 Mq, 
a ~ 0.205). The best ht based on the Hartle-Thorne geometry gives for the Gandolh 
EoS x^ - 64, a slight decrease of the mass parameter to M = 2.10Mo and a very 
slight increase of the spin parameter to a = 0.208. Therefore, we can conclude that 
the RPl-i-TPl model with the Gandolh EoS satishes the self-consistency test, as both 
precision of the ht and the estimate of the neutron star parameters are in very good 
agreement with the predictions of the Kerr approximation used in the htting procedure. 
The precision of the mass estimate is on the level of one percent. The best ht parameter 
qjcp- - 1.11 is low enough to enable the high coincidence of predictions of the Hartle- 
Thorne geometry and the Kerr approximation. 

For the RPl-i-TPl variant of the RS model and the Skyrme Gs EoS, the results of 
the htting procedure are presented in Figure|6] We again give for comparison the results 
of the best ht obtained for the Kerr geometry approximation (x^ - 61, M ~ 2.12Mo, 
a ~ 0.201); however, for this EoS the maximal allowed mass (M ~ 2.11 Mq) is slightly 
lower than the Kerr approximation estimate. The best ht based on the Hartle-Thorne 
geometry gives for the Skyrme Gs EoS x^ - 64, and a slight decrease of the mass 
parameter to M - 2.10Mo, still very close to the maximal allowed mass for this EOS, 
and a slight increase of the spin parameter to a = 0.211. Therefore, we can conclude 
that the RPl-tTPl model with the Skyrme Gs EoS satishes the self-consistency test, as 
both precision of the ht and the estimate of the neutron star parameters are in very good 
agreement with the predictions of the Kerr approximation. The precision of the mass 
estimate is on the level of one percent again. At the best ht, the parameter qja^ - 1.83 


®The fitting of the data can be done for the Hartle-Thorne geometry by considering the neutron star pa¬ 
rameters M. a, q as free parameters. However, such a fitting is extremely time consuming. We use the fitting 
tied to the EoS and the rotation frequency of the neutron star, along the curve characterized by the functions 
a(M), q(M) in the space of spacetime parameters. This is much faster procedure, being quite relevant for our 
self-consistency test. 


17 












Fig. 4: The results of the fitting procedure for the RPl+RP variant of the RS model and 
the Skyrme SV EoS.;y^ dependency on M (left), a (middle), and qja^ (right). The gray curves 
correspond to the Kerr approximation. Dashed lines correspond to values of M and a from 
Tabled] 




Fig. 5: The results of the fitting procedure for the RPl+TPl variant of the RS model and the Gan- 
dolfi EoS.;y^ dependency on M (left), a (middle) and qjcf (right). The gray curves correspond 
to the Kerr approximation. Dashed lines correspond to values of M and a from Tabled 




AC/Ms a 



Fig. 6: The results of the fitting procedure for the RPl+TPl variant of the RS model and the Gs 
EoS.^^ dependency on M (left), a (middle) and qja- (right). The gray curves correspond to 
the Kerr approximation. Dashed lines correspond to values of M and a from Table|3] 
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is still low enough to enable the high coincidence of predictions of the Hartle-Thorne 
geometry and the Kerr approximation. However, the predicted mass of the neutron star 
is very close to the maximum related to the Skyrme Gs EoS, bringing some doubts 
on the applicability of this EoS for the source 4U 1636-53; if the Skyrme Gs EoS 
is the proper one, than we could expect some strong instability of this source in near 
future. 

The results of the;if^-test realized for the external Hartle-Thorne geometry with pa¬ 
rameters governed by the acceptable EoS are summarized in Table|4] These results give 
the self-consistency test of the results obtained due to the assumption of Kerr approxi¬ 
mation of the neutron star external spacetime. We can see that only the RPl-i-TPl vari¬ 
ant of the RS model can be considered as surviving the self-consistency test. Moreover, 
the Gandolfi EoS can be considered as the most plausible one in explaining the fitting 
of observational data related to the twin HE QPOs observed in the 4U 1636-53 source. 
Notice that the best fits implied by the Hartle-Thorne geometry give in all considered 
cases value that is higher (worse fit) than in the case of the fits based on the Kerr 
approximation. Of course, we could obtain better fits by the Hartle-Thorne geometry 
for other values of the neutron star parameters. But in such a case, the fits have to 
be related to the parameters M, a, q considered as free parameters, while in our case 
the spacetime parameters were confined by the chosen EoS and the observed rotation 
frequency of the source. Therefore, we cannot exclude that in future an EoS will be 
discovered that will enable to obtain better fit to the observational data than those pre¬ 
sented in our paper. However, we can note that all EoS considered in our paper give 
fits worse than those implied by the Gandolfi EoS. 


Table 4: Results of the -test realized for the external Hartle-Thorne geometry with parameters 
governed by the acceptable EoS. 


EoS 

Models 

M [Mo] 

a 

q/a^ 


Gandolfi 

RPl-i-TPl 

2.10 

0.208 

1.77 

64 

Gs 

RPl-tTPl 

2.10 

0.211 

1.83 

64 

SV 

RPI-hRP 

2.11 

0.286 

2.99 

101 


7 Discussion 

The atoll source 4U 1636-53 seems to be one of the best test beds for both the mod¬ 
els of strong gravity phenomena and the microphysics determining equations of state 
governing the internal structure and exterior of neutron stars. This is due to simulta¬ 
neous availability of relatively good observational data of the HE QPOs occurring in 
the innermost parts of the accretion disc where the extremely strong gravity is relevant, 
enabling thus to put precise restrictions on the neutron star external spacetime parame¬ 
ters in the framework of the RS model, and the knowledge of the rotational frequency of 
the neutron star that enables a precise modeling of the internal and external structure of 
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rotating neutron stars in the framework of the Hartle-Thorne theory, for whole variety 
of the equations of state. Strong restrictions on the acceptable versions of the RS model 
can be obtained, because the precise knowledge of the rotational frequency of the neu¬ 
tron star implies a narrow evolution line for the Hartle-Thome models in the M-a 
plane that has to be adjusted to relatively precisely determined points in the M-a plane 
predicted by the acceptable variants of the RS model of the observed HF QPOs. 

The RS model can be well tested for the atoll source 4U 1636-53 since this source 
demonstrates two resonant radii in the observational data. For all relevant pairs of 
the oscillatory frequency relations of the RS model, the range of allowed values of 
th e mass and dimension less spin of the neutron star at 4U 1636-53 has been given 
Stuchlrk et alJ ( 2012 ). The most promising frequency pairs predicting the range 


in 


of the 4U 1636-53 neutron star mass and spin in accord with neutron star structure 
theory were tested by fitting the frequency relation pairs to the observational data on 
the HF QPOs separated into two parts related to the pair of frequency relations - only 
the frequency relations containing geodesic orbital a nd epicyclic frequenc ies, or some 
combinations of these frequencies, were considered ( Stuchlrk et alJ 201^ . Neverthe¬ 
less, it should be noted that the cause of the switch of the pairs of the oscillatory modes 
is not tied to the resonant phenomena related to the oscillations governed by the fre¬ 
quencies of the geodesic motion. The switch can be related, e.g., to the influence of 
the magneti c held of the neutr on star and after the switch the Alfven wave model can 
be relevant ( Zhang et al. 20061) . However, limiting the study to the resonant phenom¬ 
ena and frequencies of the geodesic origin, the number of free parameters of the model 
is restricted to the mass M and dimensionless spin a of the neutron star, as we are able 
to demonstrate that the predicted mass of the neutron star is large enough to guaran¬ 
tee with high precision independence of the geodesic frequencies on the quadrupole 
moment of the neutron star and a pplicability of the K err approximation in decribing 
the neutron star external geometry ( Urbanec et al .l2013h . Inclusion of the non-geodesic 
oscillation modes and non-resonant causes of the switch is postponed to future studies 
and could be relevant for some other sources containing neutron stars. 

The fitting procedure realized in the framework of the RS model is shown to be 
more precise by almost one order in comparison to the standard fitting b ased on the in¬ 
dividu al frequency relations that were used in pairs in the RS model ( Stuchlrk et alJ 


20141) . For example, the fitting by the standard RP model predicts the best fit along 


the mass-spin relation M{a) — Mq |^1 + 0.75{a + a^)j with rather poor maximal preci¬ 
sion of the tsst given by ~ 350 and y^/dof ~ 16; the other frequency relations 
give comparable poor precision ( Torok et alJ 2012|) . Similar results with poor preci- 
sion were obtai ned also for models with frequency relations of non-geodesic origin 
(ILin et al.ll201 ih . On the other hand, the best fit obtained for the RS model with fre¬ 
quency relation pair RPl-i-RP gives ~ 5 5 andy^/dof ~ 2.6 that is quite acceptable 
due to the character of the data distribution ( Torok et alJl2012 ). The RPl-i-TPl version 
of the RS model predicts the second best fit with precision that is given by ~ 61 and 
X^ldof ~ 2.9. 

Testing the RS model by the Hartle-Thorne theory with fixed rotation frequency 
and a variety of equations of state brings another efficient selection of tfie variants of 
the RS model. The results are illustrated by Figures [1] and and clearly demonstrate 
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that only the two variants of the RS model giving the best results of the data fitting are 
acceptable by the Hartle-Thorne models of the neutron star structure, if the rotation 
frequency /mt ~ 580Hz. The RPl+RP version of the RS model predicts mean val¬ 
ues of mass M ~ 2.20 M© and spin a ~ 0.27 and these are the data that can be met 
precisely by the Hartle-Thorne models - namely for the Skyrme equation of state SV. 
The RP1 h-TP 1 version of the RS model predicts the mean values of mass M ~ 2.12 Mq 
and spin a ~ 0.20. This mass and spin can be explained by the Hartle-Thorne model 
with the Gandolfi equation of state. It is interesting that the prediction of the Hartle- 
Thorne model based on the Skyrme equation of state Gs enters the allowed range of 
the mass and spin parameters given by the RPl -i-TPl version of the RS model, although 
it does not reach the mean value of the mass parameter, as demonstrated in Figure 
If this equation of state is the real one, the neutron star at the source 4U 1636-53 has 
to be in a state very close to instability leading to some form of collapse and dramatic 
observational phenomena. Mass and spin of the neutron star predicted by the other ver¬ 
sions of the RS model acceptable due to the data fitting are completely out of the range 
of the M-a dependencies predicted by the Hartle-Thorne model for the whole variety 
of available equations of state. 

In the special situations related to accreting neutron stars with near-maximum masses, 
the Kerr metric can be well applied in calculating the orbital and epicyclic geodesic fre¬ 
quencies, as has been done in the present paper. It should be stressed that the neutron 
star mass and spin parameters predicted by the two relevant findings of frequency pairs 
are in agreement with the assumption of near-maximum masses of the neutron stars - 
see Figure|2] For each acceptable equation of state and the observed rotation frequency 
of the neutron star, the Hartle-Thorne model has been constructed, giving thus not 
only mass M and spin a, but also the dimensionless quadrupole moment q and other 
characteristics as the equatorial radius and the compactness. The detailed results of 
the Hartle-Thorne model obtained for the three equations of state that can be in the play 
are shown in Table [3] The results clearly demonstrate that in all three cases we obtain 
a very compact neutron star, especially for the Gandolfi and Skyrme Gs equations of 
state, related to the second best fit with mass parameter M ~ 2.12 Mq and spin a = 0.2, 
having radius R ~ 3.5 M. The spin exactly predicted by the Hartle-Thorne model is 
slightly overcoming the mean value of the spin of the data fitting of the RS model, but 
it belongs to the allowed range. The parameter qla^ ~ 1.7 corresponds to the external 
Hartle-Thorne spacetime that is very close to the Kerr spacetime. 

The self-consistency test of the RPl-i-TPl variant of the RS model by the Hartle- 
Thorne geometry, related to the Gandolfi EoS, confirms this choice, since the results 
of the ^^-test are very close to those obtained due to the test by the Kerr geometry 
approximation - both for the value of;^^^ at the best fit, and the close values of the mass 
and spin spacetime parameters. Similar results are obtained by the self-consistency test 
for the Hartle-Thorne model related to the Skyrme Gs EoS. However, the test confirms 
also the conclusion that the estimated mass has to be very close to the maximum al¬ 
lowed by the EoS, lowering thus the potential relevance of this EoS for the 4U 1636-53 
neutron star. 

Eor the RPIh-RP variant of the RS model, the Hartle-Thorne model using the 
Skyrme SV EoS, related to the neutron star with mass M ~ 2.20 Mq, the spin is 
also predicted with the high precision a - 0.272, but the neutron star is not so ex- 
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tremely compact, having radius R ~ 4.1 M; the parameter qja^ ~ 2.6 is too high 
to approve the application of the Kerr geometry in description of the Hartle-Thorne 
external spacetime. In fact, the self-consistency test by the Hartle-Thorne geometry 
predicts the best fit with being too high (twice the estimate due to the fitting pro¬ 
cedure using the Kerr approximation) to imply relevance of this variant for the chosen 
EoS. For this reason the RPI-fRP variant of the RS model can be considered to be 
excluded by the self-consistence test. 

8 Conclusions 

We can conclude that there is a strong synergy effect of our approach. As expected, 
the equations of state applied in the Hartle-Thorne model of neutron stars fully exclude 
a lot of variants of the RS model that could be acceptable due to the fitting procedure to 
the HE QPO data observed in the 4U 1636-53 source. Moreover, the results of the RS 
model allow for the rotation frequency of the 4U 1636-53 neutron star /rot = 580 Hz, 
but fully exclude the possibility of /ot = 290 Hz. 

The restrictions work effectively in the opposite direction too - the results of the RS 
model put significant restrictions on the relevance of the equations of state. The crucial 
point is that the self-consistency test by fitting the observational data by the RS model 
with the orbital and epicyclic frequencies in the Hartle-Thorne geometry related to 
the acceptable EoS excludes one of the variants of the RS model predicted by the fitting 
under the Kerr approximation of the neutron star external geometry, and also the corre¬ 
sponding EoS. In fact, we have shown that the RPI-hRP variant of the RS model along 
with the Skyrme S V EoS connected to this variant are excluded by the self-consistency 
test giving high value of;^^^ at the best fit. It is interesting that this happens for the vari¬ 
ant of the RS model giving the best fit to the data of twin HE QPOs when the Kerr 
approximation of the oscillatory frequencies has been used. 

On the other hand, the RPl-i-TPl variant of the RS model related to a very stiff 
Gandolfi EoS goes successfully through the self-consistency test by the Hartle-Thorne 
geometry. In this case, the resulting best fit gives x^ - 64 that is well comparable to 
the result obtained for the Kerr approximation (x^ - 61). The resulting neutron star 
parameters (M ~ 2.10Mo, a ~ 0.208, and qja^ ~ 1.77) are also very close to those 
obtained in the Kerr approximation, demonstrating errors of one percent. Moreover, 
the Skyrme Gs EoS used in the self-consistency test by the Hartle-Thorne geometry 
gives also acceptable results, implying a possibility of the 4U 1636-53 neutron star 
being near the marginally stable state with mass M ~ 2.11 Mq. Of course, the vicinity 
of an instability of the neutron star puts some doubts on the applicability of the Skyrme 
Gs EoS. 

We can conclude that in the framework of the Hartle-Thorne theory the EoS imply 
strong restriction on the RS model of the twin HE QPOs observed in the atoll source 
4U 1636-53. In fact only the RPI-i-TPl variant of the RS model satisfies the self- 
consistency test. Moreover, it seems that there is only one EoS, namely the Gandolfi 
EoS that can be considered as a fully realistic choice in the framework of the mod¬ 
elling the twin HE QPOs. The self-consistency test also demonstrates that the Kerr 
approximation of the neutron star external geometry gives very precise estimates for 
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very compact neutron stars having sufficiently low values of the parameter qja^ < 2. 

It was shown that observations of the twin HF QPOs provide tests on equation of 
state that put limits on the gravitational mass, and the spin that is linearly related to 
the moment of inertia of the neutron star. This could provide another test of the equa¬ 
tions of state that allow for existence of neutron stars with M > 2.0 Mq. 

For the equations of state acceptable by the RS model we can determine also 
the quadrupole moment and the shape of the neutron star surface governed by the equa¬ 
torial and polar radii. These quantities have to enter other strong gravity tests of 
the 4U 1636-53 neutron star spacetime parameters predicted by the twin FIF QPOs, 
e.g., the profiled spectral lines generated at the neutron star surface or at its accretion 
disc. We believe that such tests could confirm or exclude one of the two EoS implied 
by the acceptable variant of the RS model. 

Of course, it will be very important to test the RS model of the twin HF QPOs and 
all its consequences for some other neutron star system. We have to check, if the same 
variant of the RS model, and the same EoS in the Hartle-Thorne theory of the neutron 
stars could be relevant. However, no source similar to the 4U 1636-53 neutron star sys¬ 
tem has been observed. Such sources have to demonstrate sufficiently extended range 
of the twin HE QPOs and an indication of two clusters of the observational data that 
could be related to different models of twin HP QPOs that could be switched at a res¬ 
onant point. Simultaneously, we have to know the rotation frequency of the neutron 
star. 
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Appendix. Orbital and epicyclic frequencies in Hartle- 
Thorne geometry 

Circul ar and epicyclic geodesic m otion in the Hartle-Thorne geometry has been stud¬ 
ied in Abramowicz et al. (2003); Torok et al. ( 2008b . 2015 ). Here we only present 
the expressions for the orbital (Keplerian) frequen cy and the radial and vertical 
epicyclic frequencies as given in Torok et al. ( 2008bl) . Alternative, but equivalent, ex¬ 
pressions can be found in Boshkavev et al.l (l2014l) . 

The Keplerian frequency is determined by the relations 


VK(r; M, a, q) = 
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IttGM r3/2 


1 - — + a^Ei(r) + qE2{r) 
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The radial epicyclic frequency Vr and the vertical epicyclic frequency vg are given by 
the relations 


v^(r, M, a, q) - 
Vg(r; M, a, q) - 


2nGM 


^^-[1 + aFi(r) - a^F 2 (r) - qF^ir)], (23) 
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